Evidence of two viscous relaxation processes 
in the collective dynamics of liquid lithium. 
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New inelastic X-ray scattering experiments have been performed on liquid lithium in a wide wavevec- 
tor range. With respect to the previous measurements, the instrumental resolution, improved up to 
1.5 meV, allows to accurately investigate the dynamical processes determining the observed shape 
of the the dynamic structure factor, S(Q,u>). A detailed analysis of the lineshapes shows the co- 
existence of relaxation processes with both a slow and a fast characteristic timescales, and therefore 
that pictures of the relaxation mechanisms based on a simple viscoelastic model must be abandoned. 
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The collective dynamics of liquid alkali metals exhibits 
several features which make these systems excellent can- 
didates to test different theories on the microdynamics 
of the liquid state. In simple liquids, both generalized 
kinetic theory and mode-coupling treatments predict the 
existence of two distinct relaxation processes affecting 
the dynamics of density fluctuations M. In these the- 
ories the dominant damping mechanism is provided by 
a fast process which is thought to be associated with 
the interactions between an atom and the "cage" of its 
nearest neighbors. In addition, one should also detect a 
much slower process involving the cooperative motion of 
a larger number of particles, which rearranges the local 
structure to ultimately restore equilibrium, and which 
is responsible for the critical slowing down in those sys- 
tems capable of supercooling. In simple liquids such as 
the molten alkali metals most of these predictions have 
been tested by molecular dynamics (MD) simulations [Q. 
In these numerical studies several phenomena have been 
detected, such as the persistence of well-defined density 
modes outside the strict hydrodynamic region and an in- 
crease of the sound velocity in the mesoscopic wavevec- 
tor region (the so-called "positive dispersion"). More re- 
cently, considerable attention has been devoted to check 
the existence of the afore-mentioned relaxation mecha- 
nisms [|§. 

Up to few years ago, the only experimental technique 
capable to access the relevant wavevector/frequency 
range was inelastic neutron scattering (INS). Studies 
were reported on rubidium |jj, cesium |6|, lithium |0,|), 
potassium [|[ and again rubidium |l(J . Probing wavevec- 
tors Q up to the first sharp diffraction peak, collective 
excitations outside the strict hydrodynamic regime were 
detected. The INS data have usually been analyzed with 
very simple models and, in the favorable cases (scatter- 
ing cross-section mostly coherent and speed of sound low 
enough to overcome kinematics restrictions), it was pos- 
sible to observe the positive dispersion of the sound ve- 



locity of the acoustic mode - one of the features predicted 
by the MD simulations. This situation is exemplified by 
the INS experiments of Bodensteiner et al. in liquid ce- 
sium H , where the quality of the data allowed even to 
test the lineshape predicted by the viscoelastic model, a 
sort of simplified relaxation mechanism with a single (av- 
erage) decay time. However, even in this case, the data 
accuracy did not allow to test whether the two (phys- 
ically quite different) mechanisms mentioned previously 
appear in the experimental determination of the S(Q,uj). 

The development of inelastic X-ray scattering (IXS) 
technique, paved the way for the experimental investiga- 
tions of the collective dynamics of liquids in the meso- 
scopic region. The aim of the present work is to re- 
port very high quality IXS data of the S(Q,uj) of liq- 
uid lithium, a case particularly unfavourable to INS but 
showing remarkable inelastic features [O. These data 
allow to accurately study the spectral shape, and to ob- 
tain the first strong evidence on the existence of a double- 
timescale decay mechanism underlying the dynamics of 
density fluctuations. The adoption of a spectral analysis 
scheme based on the generalized Langevin equation for- 
malism, enables us to extract quantitative information 
on the relaxation processes. 

The IXS experiment was carried out at the high res- 
olution bcamline ID 16 of the European Synchrotron 
Radiation Facility (Grenoble, F). The back- scattering 
monochromator and analyser crystals, operating at the 
(hhh) silicon reflections with h=9, 11, gave a total energy 
resolution of 3 meV for h=9 and 1.5 meV for h—ll. The 
wavevector transferred in the scattering process was se- 
lected between 1.4 nm _1 and 25 nm _1 by rotating a 7 m 
long analyzer arm in the horizontal scattering plane. The 
total Q resolution was set to 0.4 nm _1 . Energy scans, 
performed by varying the temperature of the monochro- 
mator with respect to that of the analyzer crystals, took 
about 180 min, and each spectrum at a given Q was ob- 
tained from the average of 2 to 8 scans, depending on the 



values of h and of Q. The sample (Goodfellows) had a 
nominal purity better than 0.001. The liquid lithium un- 
capped container was made of austenitic stainless steel, 
and a resistance heater was used to keep the liquid at 
475 K, i.e. slightly above the melting point at 453 K. 
Additional data were collected at 600 K. The 20 mm 
long sample, kept together by surface tension, was main- 
tained in a 10~ 6 mbar vacuum and loaded in an argon 
glove box. In the Q — E region of interest, empty vacuum 
chamber measurements gave either the flat electronic de- 
tector background of 0.6 counts/min (Q>8 ran -1 ), or, 
for Q<8 nm" 1 a small elastic line due to scattering from 
the chamber kapton windows (each 50 /im thick). This 
signal, after proper normalization, has been subtracted 
from the raw data. The scattered intensity I(Q,oj) was 
reduced to absolute units following a procedure described 
in Ref . |12f| , and the contribution from multiple scattering 
was negligible fl3|| . 

Typical spectra are reported in Fig. 1 (full dots). 
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FIG. 1. IXS spectra of liquid lithium at T = 475 

K at three wavevectors. Open Dots: experimental data. 
Dashed/Dotted lines: results of two fitting sessions (resolu- 
tion convoluted) with comparable \ 2 f° r a memory function 
model which includes thermal decay and a single-time viscous 
relaxation channel (viscoelastic model). Full line: same as 
before, except that now the viscous relaxation proceeds with 
two decay channels (microscopic and structural relaxation). 
Normalized \ 2 values are also reported, the expected value is 
1±0.1. 



The data analysis has been performed starting from a 
generalized Langevin equation: 

4> (Q, t) + wg(Q)0(Q, t) + f*M(Q, t-t')4> (Q, t')dt' = 



(here w$(Q) = 



[k B T/mS(Q)]Q 2 ) obtained adopt- 
ing a standard memory function approach for the 
density correlator <j)(Q,t)=<pQ(t)p-Q(0)>/<\pQ(t)\ > 
where p(t)={l/N) J2 Z e^'W, and <\p Q (t)\ 2 >=S(Q), 
the static structure factor. After solving this equation, 
the dynamic structure factor reads fj\f[ty as: 



S(Q,w) 
S(Q) 



■K- l ul{Q)M'{Q,u) 



■u>M"(Q,u>)f + [u>M'(Q,oj)Y 



(1) 



In Eq. (|l|), all the details of the microscopic interac- 
tions are embodied in the (complex) "memory function" 
M(Q, t), which is conveniently written as a sum of a lon- 
gitudinal "viscous" contribution Mi(Q,i) and of a term 
Mth(Q,t) arising from the coupling to thermal fluctua- 
tions. In the Q — > hydrodynamic regime, the memory 
function M(Q, £)=Ml(Q, t)+M t h(Q, t) can be written as 

M(Q,t) = M L (Q,t) + M th (Q,t) (2) 

= 2(r lL /nm)Q 2 6(t) + ( 7 - l)^(Q) e -^ 2t 

where tjl is the longitudinal viscosity, n the number den- 
sity, 7 the specific heat ratio, and Dt the thermal dif- 
fusivity. In the case of liquid metals, the inaccuracy of 
the expression (g) at finite wavevectors is particularly 
evident: due to the large thermal diffusivity of these sys- 
tems, the second term in Eq. (0) is in fact instantaneous, 
and the lineshape deduced from Eq. (|l|) is the so called 
"damped oscillator model" (DHO). This scheme is un- 
able to account for the quasi-elastic peaks evident from 
the spectra in Fig. 1, and thus must be discarded (for 
further details, see Ref. |12|). A reasonable way to im- 
prove the simple hydrodynamic, DHO-like, model (g) is 
to allow a non-instantaneous decay of the viscous term 
Ml(Q, t). On a practical basis, the simplest way of doing 
this is to assume multiple exponential decay laws for the 
time dependence of Mi,(Q,t): 



M L (Q,t) = Z? =1 Ai(Q)e 



-t/Ti(Q) 



(3) 



and to perform a refined analysis of the data by in- 
creasing N until optimum agreement is obtained. In a 
first step, we have taken N—l, which corresponds to the 
usual viscoelastic model. Then we tested the case N=2, 
which mimics the two-timescales theoretical predictions, 
an ansatz for Ml(Q, t) that was proposed to account for 
the details of S(Q, u>) in the MD simulations of Lennard- 
Jones liquids |lq] . Despite its arbitrary character, the 
exponential decay law in Eq. (||) has the advantage of 
yielding analytical Fourier transforms, thereby simplify- 
ing considerably the fitting procedure. 



The scattered intensity is proportional to the convolu- 
tion between the experimental resolution function R(co) 
and a "quantum-mechanical" dynamic structure factor 
S q (Q,ui), affected by the detailed balance condition: 



S q (Q, u) = pfiu/(l - e~^)S{Q, to) 



(4) 



The difference between the resulting function and the ex- 
perimental data has been minimized using a Levenberg- 
Marquardt procedure. ujq(Q) has been obtained from 
its definition, while the quantities 7 and Dt have been 
fixed to their experimental Q — > values, the latter being 
taken from Ref . [17| . Consequently, the only free param- 
eters in the fitting iterations were the decay times, Ti(Q), 
and the strengths of the viscous relaxations A?(Q). As a 
result, we have been able to fit our data all over the ex- 
plored Q range, and to find striking differences between 
the use of one or two relaxation mechanisms. The results 
(resolution convoluted) are shown in Fig. 1 for the data 
taken at the three reported Q values. The dashed/dotted 
and solid lines correspond respectively to the fit with one 
or two exponential functions. The fitting results provide 
definite evidences that the one decay-time viscoelastic 
model is untenable. On the contrary, a two timescale 
relaxation dynamics for the decay of Mi,(Q,i) success- 
fully accounts for the observed spectral shapes yielding 
an excellent fit (full line). The single-time viscoelastic 
fits (N=l in Eq. (g)), both initialized with free param- 
eters corresponding to the best values obtained in the 
two-times model, give results that can be read as a sort 
of "average" of the two processes, emphasizing in turn 
either the fast or the slow one. From this comparison 
one notices that the faster relaxation process accounts 
almost entirely for the width of the inelastic peak and 
for the broader part of the quasi-elastic region, while the 
slower relaxation is responsible for the narrow portion of 
the central peak. 
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FIG. 2. Ratio, A(Q), of the slow relaxation mechanism 
strength to the total strength. Full and open dots refer to 
T = 475 K and 600 K, respectively. 

Following the theoretical analysis of Ref. mM , we shall 
refer to the faster mechanism as "microscopic" (label /1) 
and to the slower one as "structural" (label a). In Fig. 2 



we report for the investigated temperatures the ratio 
A(Q) between the strengths of the slow process A„(Q) 
and the total strength A^(Q) + A^(Q) as found from the 
fitting procedure. The weight of the structural mecha- 
nism is seen to be largest at the smallest wavevectors 
(~ 20%), and to decrease with Q. This is in agreement 
with the expectation that the slower relaxation process 
affects less the dynamics of the systems as one increases 
the frequency to values sensibly larger than l/r a . 

The wavevector-dependence of the relaxation times is 
illustrated in Fig. 3. The microscopic time turns out to 
be nearly temperature independent, with values of t^Q) 
decreasing in the explored Q range from ~ 0.05 ps down 
to ~ 0.01 ps. The values of the relaxation time r a (Q) 
associated with the structural process are instead more 
scattered: r a (Q) is always of the same order as the in- 
verse of the energy resolution (~ 0.43 ps in the case of 
the (9 9 9) reflection), making rather difficult a precise 
determination. 
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FIG. 3. The times associated with fast (above) and slow 
(below) relaxation as deduced from the fits. Full and open 
symbols correspond to T — 475 K and T = 600 K, respec- 
tively. A precise determination of the slow structural decay 
time is affected by the timescale corresponding to the exper- 
imental resolution (here 3.3 meV; dashed line) 

Finally, in Fig. 4 we report the Q dependence of 
the sound speed ci(Q)=ui(Q)/Q where u>i(Q) is the 
(resolution-deconvoluted) peak position of the longi- 
tudinal current spectrum uj S(Q,oj), as obtained by 
the best- fitted values of the parameters. In the 
same figure we also report the Q-dependent isother- 



mal sound speed c (Q)—ujo(Q)/Q (expected to be the 
limiting value of Ci(Q) at the lowest Q accessible to 
IXS [Q ), as well as the " infinite- frequency" ve- 
locity c oc (QM 1 Lu 2 Q (Q) + Al(Q) + Al(Q)] l / 2 /Q as ob- 
taincd from the fit. In principle, the quantity c^Q) 
corresponds to an instantaneous, solid-like, response of 
the liquid and can be evaluated from an integral in- 
volving the effective interparticle potential and the pair 
distribution function |0,f2[- To better enlight the dif- 
ferent relaxation steps, Fig. 4 also reports the values 
of c 0O _ Q (Q) = K 2 (Q) + A2(Q)]V2/g i the » unre laxed" 

sound speed associated with the presence of the only slow 
structural mechanism. At the lowest wavevectors probed 
by IXS, the effective sound velocity ci(Q) has already 
reached c OQ - a (Q); owing to the magnitude of A(Q), its 
value is, however, only about one half of the full instan- 
taneous response Coo(Q). At larger values of Q it is ap- 
parent a "positive dispersion" of q(Q), which reaches its 
maximum at Q=6 nm _1 . 
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FIG. 4. Effective sound speed Ci(Q) (•) from the peaks 
of the deconvoluted longitudinal current spectra as deduced 
from the fits. Also reported are: the ordinary hydrody- 
namic sound speed [17] (<— arrow), the isothermal speed co(Q) 

(♦) and the instantaneous speed c^{Q) ( ) as evaluated 

theoretically in Ref. [12]. Full and open symbols refer to 
T=475 K and 600 K, respectively. From the fitting param- 
eters one may also deduce: (i) i/ojq(Q) + A^(Q)/Q ( — T — ), 
the high-frequency velocity associated with the slow process, 
and (ii) the total high-frequency velocity Crx(Q) (—A—). The 
non coincidence of the latter with the theoretical Coo(Q) can 
be ascribed to the oversimplified form of the memory function 
model. 

At even larger wavevectors (Q~10 nm -1 ), the insets 
of Fig. 3 show that the condition Wi{Q)r a u(Q)>l is at- 
tained: in this case the effective velocity ci(Q) is expected 
to approach the instantaneous value c ao (Q). Fig. 4 shows 
that this is indeed the case if Coo(Q) is evaluated theoret- 
ically from the afore-mentioned integral. On the other 
hand, the fact that the values of Coo(Q), as calculated 
from the fit parameters, are larger than the "exact" ones 
is a consequence of the oversimplified form of the two- 
exponential decay model in the short-time range relevant 



for this comparison |L2] . 

Summing up, by IXS it has been possible to detect 
all the main features of collective dynamics in a simple 
liquid metal as lithium. Specifically, our data show the 
presence of three distinct decay channels of the collective 
memory function. The first one is the "thermal" mech- 
anism: as in other molten alkali metals its relevance is 
rather small, although not entirely negligible. Much more 
important is the unambiguous evidence of two well sepa- 
rated timescales in the decay of Ml(Q, t). This provides 
a convincing experimental demonstration for the simulta- 
neous presence of both fast and slow relaxation processes 
in the collective dynamics of simple liquids. 
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